Entanglement Entropy in Random Quantum Spin-5 Chains 
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We discuss the scaling of entanglement entropy in the random singlet phase (RSP) of disordered 
quantum magnetic chains of general spin-S. Through an analysis of the general structure of the 
RSP, we show that the entanglement entropy scales logarithmically with the size of a block and 
we provide a closed expression for this scaling. This result is applicable for arbitrary quantum 
spin chains in the RSP, being dependent only on the magnitude S of the spin. Remarkably, the 
logarithmic scaling holds for the disordered chain even if the pure chain with no disorder does not 
exhibit conformal invariance, as is the case for Heisenberg integer spin chains. Our conclusions are 
supported by explicit evaluations of the entanglement entropy for random spin-1 and spin-3/2 chains 
using an asymptotically exact real-space renormalization group approach. 
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I - INTRODUCTION 

Quantum information science provides fruitful connec- 
tions between different branches of physics. In this con- 
text, the relationship between entanglement, which is a 
fundamental resource for quantum information applica- 
tions IlL and the theory of quantum critical phenom- 
ena [3, y| has been a focus of intensive research. Specif- 
ically, entanglement measures have been proposed as a 
tool to characterize quantum phase transitions (See, e.g., 
Refs. 0, H, [E B HI)- In this direction, a successful ap- 
proach has been the analysis of bipartite entanglement 
in quantum systems as measured by the von Neumann 
entropy. Given a quantum system in a pure state 
and a bipartition of the system into two blocks A and B, 
entanglement between A and B can be measured by the 
von Neumann entropy S of the reduced density matrix 
of either block, i.e., 

S = -Tr (p A log 2 p A ) = -Tr (p B log 2 p B ) , (1) 

where pA = Tr^p and ps = TiaP denote the reduced 
density matrices of blocks A and B, respectively, with 
p = \^A(tp\. By evaluating S for quantum spin systems, 
Ref. [6| numerically found that entanglement displays a 
logarithmic scaling with the size of the block in criti- 
cal (gapless) chains and saturates to a constant value in 
chains with a gap for excitations. The logarithmic scal- 
ing was proven in general for one-dimensional quantum 
models exhibiting conformal invariance [9|, , with the 
scaling governed by a universal factor, given by the cen- 
tral charge of the associated conformal field theory. In- 
deed, for a block of spins of length L in a quantum chain, 
von Neumann entropy S(L) scales as 



where c is the central charge and k is a non-universal 
constant. Recently, the behavior of entanglement en- 
tropy in critical spin chains has also been discussed in 
presence of disorder 



11, 12, 13, 14, H 



S(L) 



log 2 L + k, 



(2) 



Disorder ap- 
pears as an essential feature in a number of condensed 
matter systems, motivating a great deal of theoretical 
and experimental research (e.g., see 1(1, 17|). Moreover, 
disorder usually introduces a further effect, namely, the 
breaking of conformal symmetry in a critical model. Re- 
markably, Refael and Moore [ll| have shown that, even 
in the absence of conformal invariance, due to disorder, 
the logarithmic scaling given by Eq. holds for the 
spin- 1/2 random exchange Heisenberg antiferromagntic 
chain (REHAC) and for the Griffiths phase of the ran- 
dom transverse field Ising chain. In this case, however, 
it is governed by an effective central charge c = cln2. 
For the pure, i.e., with no disorder, antiferromagnetic 
spin-1/2 Heisenberg chain, which presents conformal in- 
variance with a central charge c = 1, any amount of 
disorder drives the ground state of this system to the so- 
called Random Singlet Phase (RSP). This is a gapless 
phase with broken conformal symmetry characterized by 
a collection of singlet pairs of spins randomly distributed 
throughout the chain. A schematic view of the RSP is 
displayed in Fig. [U The RSP appears not only for spin- 
1/2 chains, but also for higher spin disordered chains. In 
the case of integer spin chains, the RSP will usually arise 
whenever disorder is strong enough to close the Haldane 
gap. From the point of view of quantum information ap- 
plications, the RSP is potentially useful for implementing 
quantum channels with nearly perfect communication fi- 
delity, since it exhibits pairs of distant spins in highly 
entangled states [181 ]. One important open question for 
understanding the behavior of entanglement in random 
quantum chains concerns how general the logarithmic 
scaling of the von Neumann entropy is in these systems. 



FIG. 1: A schematic picture of the RSP. Spin singlets are 
composed at arbitrary distances. 



In this context, the RSP constitutes a rather interest- 
ing and universal disordered ground state, appearing in 
a number of half- integer and integer spin chains. In this 
letter, we investigate the scaling of entanglement in the 
RSP of arbitrary spin- 5 chains. We find the logarithmic 
scaling holds in general and furthermore it is associated 
with an effective central charge cr = ln(25+l) which de- 
pends only on the value of the spins in the chain. These 
conclusions apply even if the pure chain, i.e., with no dis- 
order is described by a non-conformal theory. Renormal- 
ization group calculations presented here fully support 
our results. 



II - RENORMALIZATION GROUP APPROACH 
FOR RANDOM QUANTUM SPIN CHAINS 

Renormalization group techniques provide a conve- 
nient framework for the conceptualization of the quan- 
tum phases of random spin chains. As concerns the 
RSP, it was first obtained using a perturbative real-space 
renormalization group method introduced by Ma, Das- 
gupta and Hu (MDH) [H,^ to treat the spin- 1/2 RE- 
HAC. This approach was proven to be asymptotically 
exact, which allowed for a fully characterization of the 
properties of the RSP Consider a chain of spins 

described by the Heisenberg Hamiltonian 



H — Ji~§ i ■ S 



i+l) 



(3) 



where {Si} is a set of spin-1/2 operators and {Ji} is a 
positive random variable obeying some probability dis- 
tribution P(J). The original MDH method consists in 
finding the strongest interaction between a pair of spins 
(52 and 53 in Fig. 2a) and treating the couplings of this 
pair with its neighbors ( J\ and Ji in Fig. 2a) as a per- 
turbation. The singlet formed by the spins coupled by 
the strongest bond fl is decimated away and an effective 
interaction J' is perturbatively evaluated. By iteratively 
applying this procedure, the low-energy behavior of the 
ground state is obtained as a collection of singlet pairs 
formed over arbitrary distances. This RSP is pictorially 
displayed in Fig. [TJ 

Unfortunately, when generalized to higher spins, this 
method, at least in its simplest version, revealed to be 
ineffective. The reason is that, after the elimination pro- 
cedure of the strongest bond il, the effective interaction 
J' may be greater than Q. Then, the problem becomes 
essentially non-perturbative for arbitrary distributions of 
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FIG. 2: Modified MDH renormalization procedure for spin-S 1 
chains. 



exchange interactions. For instance, considering an arbi- 
trary spin-5 REHAC, the renormalized coupling is given 
by the recursive relation [22| 



J' 



ls(s- 



1 



J1J2 
Q 



(4) 



Notice that, for 5 > 1, the renormalization factor is 
(2/3)5(5+ 1) > 1, resulting in a possible breakdown of 
perturbation theory depending on the values of J\ and 
Ji. In order to solve this problem, a gen eralization of the 
MDH method was proposed in Refs. [23, [jj] . This modi- 
fied MDH method consists in either of the following pro- 
cedures shown in Fig. [2] Taking the case of the Heisen- 
berg chain as an example, if the largest neighboring in- 
teraction to CI, say Ji, is J\ < 3/(25(5+ l))fi, then we 
eliminate the strongest coupled pair obtaining an effec- 
tive interaction between the neighbors to this pair which 
is given by Eq. (j4|) (see Fig. 2a). This new effective inter- 
action is always smaller than those eliminated. Now sup- 
pose Ji > Ji and Ji > 3/ (25(5+l))0. In this case, using 
Eq. Q would give rise to an interaction larger then those 
eliminated. In order to avoid that, we consider the trio 
of spins-5 coupled by the two strongest interactions of 
the trio, J\ and and solve it exactly (see Fig. 2b). The 
trio of spins is then substituted by one effective spin-S 
interacting with its neighbors through new renormalized 
interactions obtained by degenerate perturbation theory 
for the ground state of the trio. This method has been 
successfully applied to investigate the quantum phase di- 
agram of the spin-1 23] and spin-3/2 Q REHACs. As 
we show below, it turns out to be also essential for the 
computation of block entanglement in high spin chains. 



Ill - ENTANGLEMENT ENTROPY IN THE RSP 

Let us consider the general structure of the RSP shown 
schematically in Fig. [TJ The ground state of the sys- 
tem consists of pairs of spins coupled into singlets over 
arbitrary distances. These can be represented at ze- 
roth order by the tensor product \^p^) = \ip^ li2 ) ® 
I^m 4 > • • • ® \$i n i n+1 )> where \ip~ lk+1 ) denotes the singlet 
state (1/V2)(| Ti) — I IT)) f° r spins at arbitrary distant 
sites labelled by integer numbers i/- and ik+i- In or- 
der to evaluate entanglement, we then observe that von 
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Neumann entropy of a tensor product p ® a for density 
matrices p and a is simply S(p <£> a) = S(p) + S(a) 
Therefore, the entanglement entropy of a block of length 
L can be obtained by counting singlet pairs which cross 
the boundary of the block. This yields 



S(L)=S P (N S (L)), 



(•5) 



where S p is the entanglement entropy of a pair of spins 
and (Ns(L)} is the configurational average of the number 
of singlets connecting the two blocks. For a pair of spin-5 
particles in a singlet state, von Neumann entropy of the 
pair is maximal and given by S p — log 2 (25 + 1). Con- 
cerning (Ns(L)), it has been evaluated analytically for 
the spin- 1/2 random Hciscnbcrg antifcrromagnetic chain 
in Ref. [ll| via counting of singlets from the explicit so- 
lution of the renormalization group flow equation for the 
distribution of couplings. The value obtained for (Ns(L)) 
reads 



(N s (L)) = ^log 2 L 



(6) 



with k a non- universal constant. Although Eq. ([6]) was 
derived for the specific case of spin-1/2, it should actually 
hold for any spin-5 chain in the RSP. The reason is that 
the average number of singlet pairs should be the same for 
any ground state which is a completely random collection 
of singlets along the chain. Indeed, (Ns(L)) is a number 
which should be a general property of the structure of the 
RSP, being independent of the magnitude 5 of the spins 
composing the chain. Therefore, by inserting Eq. ^ into 
Eq. ([5]), we obtain for any random quantum spin-5 chain 
in the RSP the following relationship: 



S(L) 



cr 



log 2 L 



where 



c R = In 2 log 2 (25 + 1) = In (25 + 1) 



(7) 



(8) 



Notice that cr depends only on the value of 5, which 
is therefore the key object to set the effective central 
charge of random spin chains in the RSP. Notice that 
Eq. ([8]) does not involve any central charge of pure (non- 
disordered) models, since we do not have a central charge 
associated with integer spin chains exhibiting a Haldane 
gap. Let us turn now to the explicit evaluation of entan- 
glement entropy for the spin-1 and spin-3/2 REHACs via 
numerical implementation of the renormalization group 
equations. As it will be seen in the next sections, this 
will provide a full support to our scaling law for S(L). 



chains. The Hamiltonian for these chains will be given 
by Eq. ([3]), but with the set { Si} denoting now spin-1 
and spin-3/2 operators, respectively. Decimation of the 
chain follows the procedure previously described, with 
the renormalized couplings given by Eq. ([¥]). Then, block 
entanglement can be numerically computed in the follow- 
ing way. If a singlet is decimated and the spins compos- 
ing the singlet are in different blocks, this singlet adds 
log 2 (25 + 1) to the von Neumann entropy. On the other 
hand, in the case of a trio elimination, we have seen that 
one effective spin is returned to the chain. Then, if the 
spins composing the trio are in different blocks, the effec- 
tive spin is introduced, by a majority vote strategy, in the 
block contributing with more spins to the trio, such that, 
each block contributes with one spin for the elimination 
process. At the time the effective spin forms a singlet 
with another spin (either effective or not), entanglement 
will then be counted as in the singlet case previously 
described. Hence, notice that von Neumann entropy is 
counted only when a singlet decimation occurs. This is 
somewhat similar to the procedure adopted in Ref. 11 1 



for the transverse field Ising model, where entanglement 
was counted only for effective spins decimated by the 
transverse field. In our numerical procedure, we aver- 
aged entanglement over 1000 random coupling configu- 
rations following a power-law probability distribution of 
couplings Ji given by P(J) ~ J -0 ' 8 , for which, trio renor- 
malizations are in a negligible number in the RSP [25| . 
Numerical analysis of entanglement as a function of the 
length of the block is plotted in Fig. [31 where we consid- 
ered blocks of spins up to 10000 sites in a chain initially 
with 200000 sites. Notice that: (i) block entanglement 
scales logarithmically with the length L of the block; (ii) 
the logarithmic scaling follows our Eq. {7J , with cr = In 3 
for spin-1 and cr = In 4 for spin-3/2. The case of spin-1 
is particularly remarkable, since the scaling is governed 
by an effective central charge although the model does 
not exhibit conformal invariance even in the limit of van- 
ishing disorder. 



III. 2 - Random biquadratic spin-1 chain 

As shown by Eq. ([8]) , scaling in the RSP depends only 
on the magnitude of the spins composing the chain, being 
independent of the particular properties of the model. 
Indeed, this can be illustrated by considering the spin-1 
biquadratic chain, whose Hamiltonian is given by 



H 



(9) 



III.l - Spin-1 and spin-3/2 REHACs 

We apply now the generalized MDH method to dis- 
cuss entanglement in the spin-1 and spin-3/2 Heisenberg 



where {Si} is a set of spin-1 operators and {A;} is a 
positive random variable obeying some probability distri- 
bution P(J). Application of the generalized MDH pro- 
cedure here results only in the formation of singlets, with 
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FIG. 3: Logarithmic scaling of entanglement entropy for the 
spin-1 and spin-3/2 REHACs in the RSP. Blocks of spins up 
to 10,000 sites are considered in a chain initially with 200,000 
sites. The logarithmic fittings shown is the best curve fitting 
numerically obtained. 



FIG. 4: Logarithmic scaling of entanglement entropy for the 
random biquadratic spin-1 chain and the spin- 1/2 REHAC in 
the RSP. Blocks of spins up to 10,000 sites are considered in 
a chain initially with 200,000 sites. The logarithmic fittings 
shown is the best curve fitting numerically obtained. 



the renormalized exchange coupling reading [21 



A' 



2 AiA 2 
9 ci : 



(10) 



where Ai and A 2 are the neighbors to the strongest 
bond. Notice that trio renormalizations are completely 
absent here, since the renormalized couplings generated 
are always smaller than the energy scale CI. Similar be- 
havior occurs for the spin- 1/2 Heisenberg chain, with 
J' = (l/2)JiJ 2 /n and, therefore, the original MDH 
method can be straightforwardly applied. In this sim- 
pler situation, entanglement between a block of length L 
and the rest of the chain can be computed by exclusively 
singlet renormalizations (there is no possibility of trio 
renormalizations). Block entanglements for both the bi- 
quadratic spin-1 chain and the spin- 1/2 Heisenberg chain 
are shown in Fig. |4j As before, we considered blocks of 
spins up to 10000 sites in a chain initially with 200000 
sites, with entanglement averaged over 1000 random cou- 
pling configurations following a power-law probability 
distribution of couplings Jj given by P(J) ~ J -0-8 . For 
the spin-1/2 Heisenberg chain, the result for the scaling 
is in agreement with Ref. [llj and reproduces Eq. |5|). 
Moreover, observe that both the spin-1 REHAC in Fig. [3] 
and the random biquadratic chain exhibit the same ef- 
fective central charge cr, in complete agreement with 
Eq. ©. 



IV - CONCLUSION 

In summary, we have characterized the logarithmic 
scaling of entanglement entropy in the RSP of random 
quantum chains of general spin- S. In particular, we have 



shown that an effective central charge governing the scal- 
ing can be defined, which depends exclusively on the 
magnitude S of the spin, as given by Eq. ([5]). In or- 
der to give support to these conclusions, we evaluated 
entanglement for different spin-1 and spin-3/2 chains us- 
ing a real space renormalization group approach. These 
results are encouraging for the pursuit of a complete char- 
acterization of entanglement entropy in general random 
critical spin-S* chains. In particular, characterization of 
the complete phase diagram, i.e., for different degrees of 
disorder of spin-1 and spin-3/2 REHACS, is challenging. 
This topic is left for a future investigation. 
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Power-law probability distributions are numerically con- 
venient (but not fundamental) to make evident the scal- 
ing of entropy in the RSP, since trio renormalizations get 
rather negligible in this scenario. Indeed, our results were 
checked for a number of negative exponents, even for very 
small ones and also zero (square distribution case). 



